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Abstract 



Single particle basis functions for composite fermions are obtained from 
which many-composite fermion states confined to the lowest electronic Lan- 
dau level can be constructed in the standard manner, i.e., by building Slater 
determinants. This representation enables a Monte Carlo study of systems 
containing a large number of composite fermions, yielding new quantitative 
and qualitative information. The ground state energy and the gaps to charged 
and neutral excitations are computed for a number of fractional quantum Hall 
effect (FQHE) states, earlier off-limits to a quantitative investigation. The 
ground state energies are estimated to be accurate to ~ 0.1% and the gaps 
at the level of a few percent. It is also shown that at Landau level fillings 
smaller than or equal to 1 /9 the FQHE is unstable to a spontaneous creation 
of excitons of composite fermions. In addition, this approach provides new 
conceptual insight into the structure of the composite fermion wave functions, 
resolving in the affirmative the question of whether it is possible to motivate 
the composite fermion theory entirely within the lowest Landau level, without 
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appealing to higher Landau levels. 
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I. INTRODUCTION 



Interacting electrons confined to two dimensions and subjected to a strong magnetic 
field exhibit spectacular phenomena, e.g., the fractional quantum Hall effect (FQHE) [|l|. 
An elegant and succinct qualitative explanation of these phenomena is given in terms of new 
particles called composite fermions , which are electrons bound to an even number of vor- 
tices of the many body wave function. The reason for the simplicity of the composite fermion 
(CF) approach is that the strongly correlated liquid of electrons is mapped into a weakly 
interacting gas of composite fermions, enabling a conventional single-particle description of 
the FQHE, similar in spirit to that in Landau's theory of Fermi liquids. The composite 
fermions have the same charge, spin and statistics as electrons, but they differ from the elec- 
trons in the fundamental non-perturbative aspect that they experience a reduced effective 
magnetic field, given by 

B^^ = B- 2mp<j)o (1) 

where B is the external magnetic field, p is the electron (CF) density, and 0o = hc/e is the 
flux quantum. This is a consequence of a partial cancellation of the Aharonov-Bohm phases 
by the phases generated by the vortices on other composite fermions. The effective filling 
factor of composite fermions, u'"^ = p(j)o/\B'-''^\ is therefore given by 

„CF 

(2) 



where u = p<f)o/B is the filling factor of electrons. The minus sign corresponds to the case 
in which B'^'^ and B are antiparallel. 

The model of weakly interacting composite fermions has explained numerous facts per- 
taining to the FQHE [|l|] , some of which are as follows. The fractions appear in experiments 
in certain prominent sequences, given by 

u = , (3) 

2pn ± 1 ^ ' 
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which are called the principal sequences of fractions. These correspond in the CF theory to 
filled Landau levels (LLs) of composite fermions (z/*-^^ = n), i.e., the FQHE of electrons is 
thus simply a manifestation of the IQHE of the composite fermions. The absence of other 
fractions is related to the scarcity of FQHE in higher LLs. The experimentally determined 
gaps to charged excitations also depend on n roughly in a way expected from modeling 
them as the cyclotron energy of the composite fermions |^. An interesting limiting case is 
at z/ = 1/2, where B^^ vanishes. Halperin, Lee and Read have proposed that the composite 
fermions exist here as well and form a Fermi sea @], which has given rise to numerous 
new experiments P] that detect composite fermions by observing their cyclotron orbits 
in the vicinity of half filled Landau level and provide further nontrivial confirmation of 
the CF hypothesis. The CF approach can straightforwardly be extended to include spin 
as well 0,^, which is of interest since the Zeeman energy is quite small in GaAs. Du 
et al. have studied in detail the appearance and disappearance of the fractions as a 
function of the Zeeman energy in terms of transitions between incompressible states of 
various polarizations, originating from a competition between the cyclotron and the Zeeman 
energies of the composite fermions. Thus, the CF theory appears to give a simple and unified 
description of all the liquid states of electrons in the lowest LL, incompressible or otherwise. 

Equally striking is the success of the microscopic CF theory. The CF wave functions have 
been compared to the exact wave functions (obtained for a finite number of electrons, usually 
fewer than 10-12, by a brute force numerical diagonalization of the Hamiltonian in the lowest 
LL subspace) in a large number of studies and found to have an overlap of essentially 100% 
PJ^J^, which is particularly impressive considering that the CF wave functions have no 
adjustable parameters for the incompressible FQHE states and their low-energy excitations. 
[In contrast, in studies of atoms involving Jastrow factors, which typically depend on several 
parameters, typically only 80 to 85% of the correlation energy is obtained correctly. With 
backflow corrections and by including generalized Jastrow factors that account for three- 
body interactions, the accuracy can be increased to within a few %, but only at the cost 
of introducing many more adjustable parameters (see, for example Ref. ll^).] Therefore, 
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it should in principle be straightforward to make accurate quantitative predictions for the 
experimentally measurable quantities using the CF wave functions. 

Unfortunately, in the past, it has not been possible to compute with the CF wave func- 
tions for systems containing more than ~ 10 composite fermions due to the following tech- 
nical complication. In the simplest (unprojected) form, the CF wave function have some 
amplitude in higher electronic LLs. Even though it turns out to be surprisingly small ||11|| , 



the unprojected wave functions are not very useful for quantitative estimates, especially in 
the limit of large B. Here, a good quantitative account is given by the lowest-LL-projected 
CF wave functions. (Here and below, it will be understood that the phrase "lowest Landau 
level" refers to the lowest LL of electrons; composite fermions themselves will be allowed 
to occupy any number of CF-LLs, which all reside completely within the lowest electronic 
LL.) In our earlier work, the projection was computed numerically (exactly) by expanding 
the CF wave function into electronic basis states and throwing away the part containing 
higher LLs p[. However, the time and memory requirements increase exponentially for this 
method, limiting the projection to approximately 10 particles. It is also possible to express 
the projection operator in an integral form, but Monte Carlo with these wave functions runs 
into the technical difficulty of the fermion sign problem. 

As a result, our quantitative understanding of the FQHE in the past has been largely 
limited to what we can learn from the studies of small systems, either from exact diagonal- 
ization or using the CF wave functions (both of which produce virtually identical results). 
The size of the many-body Hilbert space is the smallest closest to half filled LL, and grows 
quickly away from u = 1/2. As a result, the maximum number of electrons for which exact 
diagonalization has been performed depends on the filling factor; at u = 1/5, 1/3, 2/5 or 
3/7, it is Nmax = 7, 8, 10, or 12 [0,115] and the largest system studied to date is 14 [Q. 
While systems containing 10-12 particles are sufficiently large for convincing ourselves of 
the validity of the concept of composite fermions and testing the accuracy of the CF wave 
functions, they are pitifully small for making reliable thermodynamic estimates for several 
quantities of interest. In fact, several of the FQHE states do not even show up in such small 



systems, which can be understood as follows. In the conventional spherical geometry, the 
n/{2n + 1) FQHE state requires at least ti^ electrons, as seen by analogy to the IQHE: 
is the minimum number of fermions needed to fill n LLs. As a result, only two data points 
are available for 3/7 (A^ = 9 and 12), while 4/9, 5/11, etc. are totally inaccessible. The 
situation is worse for the states aX v = n/ (4ri±l) - here, it is difficult to study even 2/9. Due 
to an exponential increase with the number of electrons in the computer time and memory 
requirements, it is safe to assume that at best marginal improvements will be possible in the 
future exact diagonalization studies, which will therefore not tell us much more than what 
they already have. Consquently, our best hope for a better quantitative description of the 
FQHE lies in developing new techniques to work with the CP wave functions. 

Progress was made in studying large CP systems for the low-energy excitations of the 
u = l/{2p + 1) states. Here the unprojected CP wave functions have the special feature 
that they have at most one electron in the second LL, and the projection operator can be 
written simply as K — ^hujc, where K is the kinetic energy operator. This was exploited by 
Bonesteel to study large CP systems and to compute the gap to charged excitation. The 



CP prediction for the gap was slightly lower than the one obtained from the Laughlin's trial 
wave functions. The full dispersion of the neutral excitation was determined by Kamilla, 
Wu and Jain using similar techniques ITSf. This approach, however, is not applicable to 



general fractions, for which the number of electrons in higher LLs in the unprojected CP 
wave function is a finite fraction of the total number of electrons, growing with the total 
number of electrons. 

It is clearly desirable to find ways of computing with the CP wave functions for large 
systems at arbitrary filling factors. The study of n/{2pn + 1) states at large n is interesting 
not only in its own right but may also shed further light on the compressible 1 / (2p) state. In 
this paper, we report on a method which enables computations for rather large CP systems. 
We give below results for as many as 50 composite fermions; treatment of bigger systems 
should be possible in the future. This constitutes a significant step in our ability to achieve 
detailed quantitative predictions for the PQHE. Short reports on parts of this paper have 



been published elsewhere [0,0 • 



As a first application of this method, we have computed the ground state energies and 
the gaps to charged and neutral excitation for several FQHE states of interest. There are 



several parameters in actual experiments, namely the finite quantum well width |18], Landau 



level mixing finite Zeeman energy, and disorder, which influence these energies. Our 
aim in this work will be to obtain a quantitative description in an idealized situation where 
these effects become unimportant, namely the limit in which the transverse width of the 
quantum well, disorder, and {e^ / d) / {huo^ all vanish. The last has the effect of supressing 
any LL mixing or spin reversal. Besides being theoretically convenient and clean, this limit 
is sufficient for obtaining zeroth order estimates as well as trends. Furthermore, previous 
work tells us how to include these effects and how much of a quantitative correction they 
make. For example, the gap to charged excitations is reduced approximately by a factor of 
2 when corrections due to finite well width and LL mixing are taken into account for typical 



sample parameters |20[. In fact, all of the effects neglected in the ideal limit are found to 
reduce the gaps, so our estimates will provide an upper bound for them. 

Another promising approach toward a quantitative description of the FQHE is the Chern- 



Simons field theoretic approach pT|J^. It is best suited for a description of the long distance 
physics, an understanding of which is most crucial in the case of the Fermi sea of the compos- 
ite fermions, where the correlations are expected to have power law decays. A computation 
of the short distance properties like the gaps has not yet been possible in this approach, 
except asymptotically close to z/ = 1/2 |2^. In the most widely used framework [Q, the CS 



field theory is treated as an effective long-distance theory of composite fermions, with the 
effective mass of composite fermions treated as an input parameter, to be fixed either by 
comparison with experiment or by other microscopic means. 

In addition to providing a new calculational tool, the new lowest LL (LLL) represen- 
tation of the CF wave functions also suggests how the CF theory can be motivated quite 
plausibly within the lowest Landau level, without appealing to higher LLs at all. The com- 
posite fermions, their effective magnetic field, their quasi-LL structure, and their effective 



cyclotron energy are all seen to originate in this approach as a direct consequence of repulsive 
interelectron interaction. The single particle wave functions of composite fermions in the 
lowest (electronic) LL are shown to possess the appealing property of having zeros tightly 
bound to electrons, aside from a few localized defects. This shows that while the analogy 
to higher LLs is a convenient, and intuitively powerful route to the CF theory, it is by no 
means necessary. 

The organization of the rest of this chapter is as follows. Section II gives the explicit form 
of the CF wave functions in the lowest Landau level. The essential result is displayed in the 
beginning of II. C in a form ready to be used by the interested reader. Section III presents 
results of our Monte Carlo calculations for the ground state energies and the gaps. In 
particular, it is found that the gaps to charged excitations are in reasonable agreement with 
the formula suggested by Halperin et al. 0] , and the study of the neutral excitations uncovers 
an instability of FQHE for v < 1/9. In Section IV, we show how the composite fermion wave 
functions can be motivated within the lowest LL, and how concepts like effective magnetic 
field, effective filling factor, CF-LL index, and effective cyclotron energy appear here. This 
work is formulated in the disk geometry, which also gives us an excuse to compare the CF 
wave functions with the exact wave functions in the disk geometry. The paper is concluded 
in Section V. 



II. COMPOSITE FERMION WAVE FUNCTIONS IN LOWEST LANDAU LEVEL 

We start by considering the spherical geometry in which N electrons move on the two- 
dimensional surface of a sphere [^,0] under the influence of a radial magnetic field B. By 
virtue of having no boundaries, this geometry is convenient for an investigation of the bulk 
properties of the CF state. The disk geometry will be considered in a later section. 



A. Single electron states 

The single particle eigenstates are the monopole harmonics, Yq^n,m P^ : 
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where x = cos 9, the total flux through the sphere is equal to 2g0o (2^7 is an integer), 
n = 0, 1, 2, ... is the LL index, m = —q — n,—q — n+l,...,q + n labels the degenerate states 
in the nth LL, and I = q + n. 9 and are the usual angular coordinates of the spherical 
geometry, and P"'^ is the Jacobi polynomial. The normalization coefficient is given by 

_ / (2g + 2n + 1) (g + n - m)\{q + n + m)\ V'^ 
"""^ " 1^ n!(2g + n)! ) ' ^ ^ 

Substituting the explicit form for the Jacobi polynomial, Y can be rewritten asQ 



v*vr-'(u*u) 



(6) 

where fl represents the angular coordinates 9 and of the electron, and 

u = cos(0/2) exp(-i0/2) (7) 

w = sin(^/2)exp(i0/2) . (8) 



B. Composite fermion states — unprojected 

Now we consider interacting electrons at flux Q. The composite fermion theory 
postulates that the strongly correlated liquid of interacting electrons at Q is equivalent to a 
weakly interacting gas of composite fermions at g = Q—p{N — l). (We omit the conventional 



superscript CF or * on g to simplify the formulas.) The unprojected wave function for 
the CF state at q, <|)'^'^~"p, is related to that of non-interacting electrons at q, $, as 



^Throughout this work, the binomial coefficient is to be set equal to zero if either /3 > 7 or 
/? < 0. 
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CF-up 



(9) 



the Jastrow factor J is given by 



J = = W_{ujVk - VjUkf^ exp[ip{(f)j + (j)k)] 

j<k 

where $i is the wave function of the lowest filled LL. 



(10) 



C. Lowest LL projection 

Since we are interested in the limit of i? ^ oo, we will be using the LLL projections of 
(^CF-up^ denoted by <l>^^: 



<^CF ^ p^^CF-up 



(11) 



where V is the LLL projection operator. We begin by quoting the final result |T6[ in a form 
that can be used directly by someone not interested in the details of its derivation, which 
fill the rest of the subsection. 
Final Result: 

$, the wave function of electrons at q, is a linear superposition of Slater determinant 
basis states made of Yq^n,m- The corresponding wave function of composite fermions at q 
(which gives the wave function of interacting electrons at Q) is same as $, except with Yq^n,m 
replaced by Y^^„^, given by 



q+n—m 



y2Q -|- 1)! q+m q—m 
7/ ■ H- 

(2Q+n + l)! ' ' 



7^0 \s)\q + n-m- s) ^ ^ 



d 



d 



^ duj J \ dvj 



where 



(12) 



Jj = ]li'>^j^k - VjUk) exp[-(0j + (pk) 
k ^ 



(13) 



Here the prime denotes the condition k ^ j. 
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This implies that Y^^,^{Qj) can be interpreted as the single particle basis functions of 
composite fermions, using which many-CF states can be constructed in the usual manner, 
in complete analogy to the non-interacting electron states. 

In order to deal with the derivatives in our Monte Carlo calculations, we find it convenient 
to write them as follows: 



where 



'JL 



'JL 

dVn 



J 3 



(14) 



d 



Vk 



+ 



d 



UjVk - VjUk duj 



(15) 



d 



+ 



d 



(16) 



For a given n, explicit analytical form of the derivatives is used in the evaluation of the 
wave function. The advantage of this form is that the J/s can be factored out of the Slater 
determinants to give back the usual Jastrow factor 



(17) 



Note that the phase factors exp[|(0j + 0^)] do not play any role in the calculation of the 
energies. 

Derivation: 

To derive this result, we consider a case when $ is a single Slater determinant, given by 

yi(Qi) ^1(^2) . . 



$ = Det[Yi{nj)] 



(18) 



so the unprojected CF wave function is given by 
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^CF-up ^ j]jet[Y,{nj)] . (19) 

A generalization to situations in which $ is a hnear superposition of many Slater determi- 
nants is straightforward. We need to show that ^""^ is obtained by replacing F's in Eq. (|18D 
by F^^'s. 

The Jastrow factor can be incorporated into the Slater determinant to give 

^CF^up ^ Det[Y,{nj)Jj] . (20) 

We now project this on to the lowest LL by projecting each element onto the lowest LL. I.e., 
we write 

= Det[VY,{Qj)jP] . (21) 

Since the product of two LLL wave functions is also in the lowest LL, is guaranteed to 
be in the lowest LL. It remains now to show that 

VY,,^,U^,)J^ = Y,IU^,) . (22) 

In order to evaluate the projection we first show that there exists an operator „ m 
satisfying the property that 

'T^Yq^n,mYqi fl^ra' ~ ^ q^n,rri^g' ,0,m' i (2^) 

where yg',o,m' ~ e^'^''^^u'j vj is a LLL wave function for monopole strength q' . For the 

q' 

present purposes, it is important that „ ^ be independent of m' although it may depend 
on q'; this is because the single particle wave functions of the jth electron appearing in Jj 
have different m' but the same q' = {N — l)/2. To this end, we multiply one of the terms on 
the right hand side of Eq. (||) by the LLL wave function y^'^cm' and write (with Q = q + q', 
M = m + m', and the subscript j suppressed): 

eiO0(^*^)n-s(y*y)^^Q+A/^Q-M _ age^^'^M'3+*^t;^-*^ + higher LL states. (24) 
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For \M\ > Q, must vanish, since \M\ < Q in the lowest LL. Let us first consider the 
case \M\ < Q. Multiphng both sides by e~*'3'^u*'3+^^i;*'3~*^ and integrating over the angular 
coordinates gives 

^ {Q + M + n-s)\{Q- M + s)\{2Q + 1)! 

iQ + M)\{Q-M)\{2Q + n + l)\ ' ^ ^ 

This shows that, apart from an m'-independent multiplicative constant {2Q + 1)1/ {2Q + 
n + 1)!, the LLL projection of the left hand side of Eq. (|2^) can be accomplished by first 
bringing all u* and v* to the left and then making the replacement 

^* ^ ^ ^* ^ ^ ('26) 

du ' dv ' 

While true in general for \M\ < Q, this prescription can be shown to produce the correct 
result (i.e., zero) even for \M\ > Q for the LLL projection of states of the form Vq,n,m^g',o,m', 
as proven in the Appendix. Thus, in Eq. (pSj), 



_ {2Q + 1)! ,_^.q+n-miq^j 



n—s 



A delightful simplification occurs when one brings all the derivatives to the right in 
Eq. using 

and a similar equation for the derivative with respect to u (with the summation index a'). 
The sum over s in Eq. ( p7|) then takes the form 

n—a ^ ^ /\ n— Q— a' ^ ^ /\ 

"""7 = i: (-1)"'-" 7 • 

s=a' \ s a J \ s / 

which is equal to (— 1)"'(1 — l)""""""' and vanishes unless n = a+a' . The only term satisfying 
this condition is one with a = n — s and a' = s. Consequently, the derivatives in Eq. ( P7[ ) 
can be moved to the extreme right and act only on the following LLL wave function. This 
gives Eq. (|2^), completing our proof. 
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Before ending this subsection, we note that with the help of the identity u*u + v*v = 1, 
the Eq. (|^) can be cast into different equivalent forms; for example, the sum over s can 
be expressed as a power series entirely in v*v or u*u. While the projection can be carried 
out with any form, the subsequent formulas are simplest expression in Eq. is used. 
Otherwise, the simple replacement of u* or v* by the corresponding derivatives does not 
give the projection and neither are the cancellations indicated above obtained. 



III. LOW ENERGY STATES 

The most important goal of theory is to provide an understanding of the low-lying states 
of the system in question. We will concentrate here on the ground and the low-energy states 
at the special filling factors u = n/ {2pn + 1). The ground state energy of the incompressible 
states is of interest, even though it cannot be measured directly, because it serves as an 
important criterion when comparing different theories of the FQHE, or in investigating the 
question of the stability of a FQHE state relative to another state with a different symmetry, 
e.g., the Wigner crystal p5| , |26[] . Two kinds of excitations are relevant for experiments: the 
charged and the neutral excitations. The gap to the charged excitations appears in transport 
experiments as the activation energy, while the dispersion of the neutral excitations, in 
particular the gap at the minima or maxima, can be measured by inelastic light or 
phonon scattering The smallest energy required to create an excitation will also play 
an important role in the low-temperature thermodynamics of the FQHE. 

Let us first summarize the status of our quantitative understanding of the FQHE prior 
to the the CF approach, (i) Ground state energy: The energies of the 1/3 and 2/5 states 
have been determined reasonably accurately from an extrapolation of exact diagonalization 



results |T3|J29|] . A good ansatz exists for the ground state at u = l/{2p + 1) in the form 
of the Laughlin wave functions. The energy of the Laughlin state is known accurately 
from Monte Carlo study |^ of large systems, and is in good agreement with that obtained 
from the exact diagonalization results, (ii) Charge gap: Almost the same story holds here. 
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Reasonable estimates are available for the gaps of 1/3, 2/5 and, to some extent, for 3/7 from 



exact diagonalization studies |2^. The gap of the 1/3 state has also been calculated using 
the Laughlin trial wave functions, (iii) Neutral excitations: A single mode approximation 
(SMA) [^] was found to work reasonably well for 1/3, but less well for other Laughlin 
states, as we shall see, and not well for other FQHE states. The exact diagonalization 
studies exhibit strong finite size effects for the n/ {2n + 1) states with n ^ 1 and fail to give 
a reliable account for the dispersion of the neutral excitations. 

In the CF theory, the FQHE ground state at u = n/{2pn + 1) corresponds to n filled 
CF-LLs. The excitations are obtained by removing one composite fermion from the nth CF- 
LL, leaving behind what we will call a "quasihole", and placing it in the {n + l)st CF-LL, 
creating a "quasiparticle" . This excitation will be called a (CF-)exciton. The orbital angular 
momentum L of the state is related to the wave vector k of the planar geometry as L = kR 



where R = \fQ is the radius of the sphere, in units of the magnetic length (/ = ^UcjeB). 
The configuration with the quasiparticle and quasihole on opposite poles will be identified 
with the charged excitation, as the distance between the quasiparticle and the quasihole is 
the maximum; this corresponds to the largest L state of the exciton. 

The transport gaps at z/ = 1/(2]? + 1) have been computed using the CF theory |T^,^. 
It was found that the CF theory gives marginally smaller gaps than those computed using 
Laughlin's trial wave functions (the difference is more significant for short range interac- 
tions than for the long range Coulomb interaction). We note here that here the CF wave 
functions for the ground state and the quasihole are identical to those of Laughlin's, but the 
quasiparticle wave functions are different, resulting in unequal values of the gaps. 

A. Calculation Procedure 

The interaction energy per particle is 
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Rij = 2i?| cos(^i/2) sin(^j/2) - cos(^j/2) sm(^i/2)e^('^'-'^^)| , (31) 

where Rij is the arc distance between two electrons, and R = ^/Q I is the radius of the sphere. 
It will be evaluated by variational Monte Carlo techniques, with the absolute value squared 
of the Slater determinantal wave function of composite fermions used as the probability 
weight function for the Metropolis acceptance sampling. The methods are standard, but a 
few points should be mentioned. 

It takes OiN^) operations to evaluate a Slater determinant. For a free Fermion problem, 
if a single particle is moved, only one row (or column) of the density matrix is altered, 
and it takes only operations to evaluate the new determinant in terms of the old inverse 
matrix and (9(A^^) operations to update the inverse matrix if the move is accepted p4 



However, the CF Slater determinant describes a strongly correlated electron system and 
hence is explicitly non-local, i.e. when a single electron is moved all CF matrix elements are 
altered, since the single particle wave function of a given composite fermion depends on the 
positions of all the particles. Therefore, an update of the trial N x N density matrix of the 
CF system proceeds through 0{N^) operations at every Metropolis step. 

The time T involved in our variational Monte Carlo calculations for the coulomb inter- 
action energies of the N - electron system in the ^2p2±i\ -^QH state scales as: 



T oc Nh'^n^ 

the cost of evaluating the CF determinant goes as A^^, the cost of projection as n^, and 
with e = , where e is the desired simulation error and Nmr is the total number of 



' mc 



Monte Carlo steps. (The Monte Carlo statistical error scales as (t^/ y Nmc/^c, where 
is the standard deviation in the wave function and C,c is the correlation time, which is a 
property of the wave function being sampled.) Typically, we have performed calculations 
on systems containing upto 50 particles. For example, a calculation of the ground state 
interaction energy on the z/ = 2/5 CF state for 30 particles with about 20 million Monte 
Carlo steps takes ~ 30 CPU hours on our work-station (DEC Model 250/4, CPU 21064A, 
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Mhz 266). A similar calculation for 30 particles in the z/ = 3/7 CF state takes ~ 100 CPU 
hours on the same machine. 

For the ground states and the charged excitations a.t u = n/{2pn + 1) the CF wave 
functions are single Slater determinants. For the neutral exciton, the wave function is a 
linear superposition of ~ N/n determinants, making the computations significantly more 
time consuming. In addition, the energy of the exciton has to be evaluated with extremely 
high accuracy in order to obtain reasonable values for the energy differences; up to 4 x 10^ 
Monte Carlo steps were carried out in our computations of the exciton energy. A good 
accuracy is especially important at small u, where the exciton energy becomes rather small; 
we will be dealing with energy differences as small as 0.001 e^/e/, which is approximately 
two orders of magnitude smaller than the transport gap of the 1/3 state. 

We note that when Coulomb energies are calculated for finite systems with an intent 



to obtain thermodynamic estimates, usually an infinite Ewald summation |3J] has to be 
performed that accounts for the interaction of the charges with their image charges with 
respect to the finite boundary imposed in the calculation; no such sum is necessary in the 
spherical geometry due to absence of boundaries. 



B. Accuracy of CF wave functions 

Before we proceed to compute the energies, we ask how reliable the CF wave functions 
are. To this end, we compare the energies of the CF wave functions to the exact eigenenergies 
obtained numerically, both for the ground state and the charged excitation, shown in Table 
I. The exact energies were obtained by a numerical diagonalization of the Hamiltonian in the 
lowest LL subspace, assuming an ideal two-dimensional system. Further details can be found 
in the literature (see, for example, Wu and Jain All energies quoted in Table I and the 
rest of the paper include the interaction with the positively charged neutralizing background, 
which amounts to an addition of —N^/2y/Q, the interaction energy of the N electrons with 
a positive charge +Ne at the center of the sphere. The accuracy of the CF theory is 
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impressive; the predicted energies agree with the exact energies to better than 0.05% for the 
systems studied, which represent the biggest systems for which exact diagonahzation has 
been performed at 2/5 and 3/7. We note that these system sizes are aheady large compared 
to the relevant length in the problem, namely the magnetic length (or the effective magnetic 
length), which, coupled with the fact that the energies are determined largely by the short- 
distance behavior, suggests that the error in the thermodynamic limit ought not to be much 
larger. In other words, extrapolation of the CF energies to large will produce virtually 
the same results as extrapolation of the exact energies. 

Since the gap is an 0(1) energy obtained as the difference between two large 0{N) 
energies, its accuracy is expected to be less than that of either the ground state or the 
excited state energy. A comparison of the gaps predicted by the CF theory with with the 
corresponding exact gaps (Table I) shows that the discrepancy is of the order of a few 
percent. This level of accuracy is possible only because the ground and the excited state 
energies are extremely accurate. The full dispersions of the neutral excitations obtained in 
the CF theory have a similar level of accuracy [P^. 

It is obvious that the ground state wave functions considered here contain no adjustable 
parameters as they are related to the unique filled LL wave functions. Though less obvious, 
this is also true for the single exciton state as well (at arbitrary wave vector); they are 
also completely determined by symmetry in the restricted Hilbert space of the CF wave 
functions. Also, since the CF wave functions constructed here are strictly within the lowest 
electronic LL, their energies provide strict variational bounds in the limit of zero thickness, 
no disorder, and large magnetic field. 

Another point that needs to be emphasized is that although the CF wave functions are 
analogous to, and indeed motivated by the wave functions of non-interacting electrons at q, 
they describe fully interacting electrons at Q. In other words, the CF wave functions provide 
an account of the residual interaction between composite fermions themselves. In this sense, 
the wave functions go beyond the simple mean-field picture of non-interacting composite 
fermions. We will see below that the residual interaction between composite fermions is 
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responsible for the destabilization of the FQHE for z/ < 1/9 ||17|| . 



C. FQHE Ground states 

Having estimated the level of confidence in the CF theory, we proceed to study larger 
systems. The energies of the CF ground states for several values of at 2/5, 3/7, 4/9, 
5/11, 2/9, 3/13 and 4/17 are given in Table III. While estimating the N ^ oo limit, we find 
it useful to correct these energies for the finite size deviation of the electron density from 



its thermodynamic value, p, by multiplying the energies by a factor y p/ Pn = ySQzz/A^, 
following Ref. [0. This substantially reduces N dependence and thereby allows for better 
extrapolation to the thermodynamic limit, shown in Fig. |l|. The thermodynamic estimates 
are obtained using a chi-square fitting of the energy as a polynomial in 1/A^ that biases 
the points by their error bars; the number quoted in Table III is an average over estimates 
obtained through linear, quadratic, and cubic extrapolation. 

Fig. shows the ground state energy of the CF liquid state along with the best estimate 
for the energy of the Wigner crystal It illustrates nicely the fundamental reason for 



why the Wigner crystal state is not observed as soon as the kinetic energy is suppressed 
by forcing all electrons to occupy the lowest LL: it is preempted by the CF liquid which 
has substantially lower energy. We note here that while the energies plotted in Fig. |^ 
are variational estimates, they are believed to be extremely accurate approximations to the 
actual energies, so, the conclusion drawn from this figure is reliable except in the crossover 
region near u = 1/7, where the two energies are very close. We will have more to say in this 
regard later; here we would like to stress that our estimates seem to rule out Wigner crystal 
in the vicinity of u = 1/2. 

Due to particle-hole symmetry in the lowest LL, the ground state energy per particle of 
the state at 

u' = l ^ (32) 

2pn + l ^ ^ 

is related to the ground state energy per particle of the state at ly — 2pn-\-l 
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u[E,{u) - E,il)] = (1 - u)[E,il E,{1)], 



(33) 



where Eg{l) = —yir/S, is the Hartree-Fock energy of a filled Landau level, and all energies 
are in units of e^/e/ ||36|| . 

D. Transport Gaps 

Next we come to the energy of the charged excitation (the CF exciton with the largest 
L). For z/ = 1/ {2p+ 1), the gap computed with the present projection technique is identical 
to that computed with the full projection, investigated earlier by Bonesteel Jl^, and we will 



simply use his results. The gaps for several values of at 2/5, 3/7, 4/9, 5/11, 2/9, 3/13 and 
4/17 are given in Table III. We subtract from the gap — (2p + l)^^/(2y^), the interaction 
between pointlike particles of charges e/{2p+ 1) and —e/{2p+ 1) placed at the two poles, to 
correct for the interaction between the quasiparticles; these are indicated as the corrected 
gap in Table III (no correction is made for the density). Fig. |^ shows typical dependence of 
the gap on the number of particles A^. The N ^ oo estimates of the gaps are given in Table 
III. We note that the gaps of two states related by particle-hole symmetry are the same, 
when measured in units of e^/e/. 

Apart from the intrinsic error of the CF wave functions, believed to be quite small, there 
are two sources of uncertainty in our estimates, both of which can be further reduced, if so 
needed. One is, of course, the statistical error in our Monte Carlo, and the other comes from 
the non-monotonous, oscillatory part of the interaction between the quasiparticle and the 
quasihole as a function of their distance, originating from the fact that the quasiparticles 
are not point-like but extended in space and their density exhibits oscillations as a function 
of the distance from the center, analogous to the LL wave functions. As a result, the 
interaction energy becomes equal to the correction term used above only when the distance 
is very large compared to the size of the quasiparticles, which scales as ~ v^2n times the 
effective magnetic length, requiring a study of bigger and bigger systems as we go to larger 
n along the njifln + 1) sequence. For 4/9 or 5/11, such fiuctuations can be seen even in 
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systems with as many as 40 fermions. A reduction of this error will require investigations 
of larger systems. 

In a mean-field picture, it is natural to equate the transport gap to the effective cyclotron 
energy of composite fermions. This would give 

eB^^ eB 



However, since the gap in the LLL constrained theory must be proportional to V-B, the 
effective mass of composite fermions cannot be a constant but rather must be mP^ ~ \fB. 
This kind of reasoning led Halperin et al. [|] to conjecture that the gaps of the n/{2pn + 1) 
states are given by 

C 

= (2pn ' ^^^^ 

where it was estimated from the then existing results that the constant C ~ 0.31. As the 
Fig. ^ shows, the overall behavior of the gaps computed here is in good agreement with 
Eq. ( |55D for the p = 1 sequence. When the gaps are plotted in constant units (say, Kelvin) as 
a function of (2n + 1)"^, the best straight line fit has a negative intercept. A fit as a function 
of 1/n is also not as satisfactory. Thus, our calculations provide a microscopic justification 
for Eq. (|35|). For the p = 2 FQHE states, the fit is less satisfactory. The analysis of Stern 
and Halperin |^ predicts logarithmic corrections to the formula in Eq. (|35|) , which however 
are small except in the close vicinity of the half filled LL. At the level of the accuracy of our 
study, we do not expect to be able to distinguish such small deviations. 

At the moment, it is not possible to compare the values of the gaps obtained here with 
experiments, due to the neglected effects of finite width of the quantum well Landau level 
mixing and disorder. Nevertheless, the trend predicted by the CF theory, as summarized in 
Eq. (|35|), is in good agreement with the experimental behavior. The experimental data 
are reasonably well described by 
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which suggests that the finite width and LL mixing might, to zeroth order, only alter the 
value of C, and disorder might enter in the form of a constant shift of the gaps, which was 
interpreted by Du et al. as the disorder induced broadening of the CF-LLs. 

Our work can be generalized to include some of the effects neglected in the system 
considered above. The finite well width softens the Coulomb interaction at short distances, 
which is straightforward to incorporate into our Monte Carlo and will be the subject of a 
future work. Landau level mixing is harder to deal with. It has been investigated by using 
Jastrow factors that are not analytic in Zj and thereby contain some LL mixing. Fixed phase 
Monte Carlo ||3^ or a variational approach considering a linear combination of <|)*-^'^~"p and 
13811 will be useful in this context. 



E. Neutral excitons 



Next we come to neutral excitations |T5|. Consider a state with a hole in the rrih state 
of the topmost occupied CF-LL and a composite fermion in the me state of the lowest 
unoccupied CF-LL. With no loss of generality, we assume = for the excited state, 

i.e. m/j = me, and denote the corresponding Slater determinant basis state by |m/j >. Then, 
the exciton state with a definite total angular momenta L is a linear superposition of |m/j > 
with different m/^, given by 

^exciton = J2<nt + q, -ruh] nt + l + q, mh\L, > [m^ > , (37) 

where rit is the LL index of the topmost occupied CF LL, and nih = —q — rit, ...,q + nt. Note 
that the coherence factors are the same as those for the IQHE exciton. 

Figs. ||-^ show the energy of the exciton, A^^', at z/ = 1/5, 1/7, 1/9 and 2/5 (measured 
relative to the ground state energy) as a function of the wave vector k = Lj^/Q. We have 
corrected both the ground and the exciton energies for finite size deviation of the density, as 
discussed earlier. The fact that the results for systems with different numbers of electrons fall 
on a single line shows that remaining finite size effects are small. For strictly non-interacting 
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composite fermions, one would expect no dispersion of the exciton energy; the dispersion is 
a measure of the interaction between the quasiparticle and the quasihole as their relative 
separation is varied. The oscillatory structure arises because of the internal oscillations in 
the density of the quasiparticles. 

The most striking result is that the energy of the CF exciton at z/ = 1/9 falls below that of 
the ground state at approximately kl ^ 0.85. The FQHE is thus explicitly demonstrated to 
be unstable to a spontaneous creation of CF excitons. Since the size of the minimum energy 
exciton {kl^/u ^ 7.51) is approximately equal to the interparticle separation, a large number 
of excitons can be created before their interaction may limit their further production. This 
will clearly lead to a collapse of the FQHE, signaling the destruction of composite fermions 
through a vortex unbinding transition. [We note that there is another mechanism for a 
transition to the Wigner crystal state: for u < 1/72 the Laughlin wave function itself 
describes a Wigner crystal |3^. Here the composite fermions are not destroyed, i.e., the 
FQHE liquid of composite fermions freezes into a Wigner crystal of composite fermions.] 
Furthermore, the fact that the wave vector of the instability is close to the reciprocal lattice 
vector of the WC at z/ = 1/9, kwcl = {^t^v j ^pifl'^ ^ 0.898 also gives a hint to the 
nature of the true ground state. Similar result will undoubtedly be obtained for v < 1/9, 
ruling out FQHE for z/ < 1/9. 



Earlier studies |^5[ have compared the energy of the Laughlin state with that of the 
Wigner crystal state, and predicted a transition at z/~^ = 6.5 ±0.5 (see. Fig. ^for example). 
In particular, it has been well established that the incompressible FQHE state is not the 
ground state at 1/9. However, this is the first time that an instability of a FQHE state has 
been found, despite a number of previous studies investigating the neutral excitations pl|,p!2|. 



(It should be noted, however, that the field theoretical approaches of may in principle 
be capable of ruling out FQHE by finding either instabilities or other lower energy saddle 
point solutions.) It is remarkable that this can be discovered entirely within the tightly 
constrained, zero-parameter framework of the composite fermion theory. We note that this 
instability is also an explicit demonstration of the breakdown of the mean-field description 
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of the CF state, as a result of level crossings in going from the mean field description to the 
true state. 

At V = 1/7, the minimum energy of the exciton, A^j„, becomes rather small but remains 
positive. This does not rule out a WC ground state here, since a first order phase transition 
to a WC state may take place before A^-^ hits zero. Nonetheless, our result indicates 
that the question of which state is the ground state at z/ = 1/7 deserves a more careful 
analysis. Comparison with Laughlin wave function favors a WC ground state here. 



but the Laughlin wave function becomes progressively worse as one goes to smaller filling 
factors and may not be sufficiently accurate to lead to a definitive conclusion at u = 1/7. 



To elaborate, in the numerical results of Ref. |13], the energy of the Laughlin wave function 



at z/ = 1/5 is found to deviate by ~ 0.15% from the exact energy already for 7 electrons. 
The error in the thermodynamic limit will be more and for u = 1/7 still greater. While 
the error is quite small in absolute terms, it may be significant while comparing with the 
Wigner crystal state. Our accurate estimate for the energy of the Laughlin 1/7 state is 
-0.280944(6) e'^/el while the best estimate for the WC energy is -0.2816(5) e^/el fS^. The 
difference is ~ 0.2%, which is likely to be small compared to the amount by which the 
Laughlin state overestimates the energy of the FQHE liquid at 1/7 in the thermodynamic 
limit. On the other hand, the energy of the WC state can also surely be improved. Thus, 
further theoretical work will be needed to provide a conclusive answer to which state has 
lower energy. It is likely that the issue will eventually be settled by experiments. In this 
context, it is worth pointing out that in experiments in semiconductor heterj unctions a 
significant reduction of the longitudinal resistance (15-20%) at z/ = 1/7 has been reported, 
providing evidence of FQHE at 1/7 (not conclusive though, as no plateau has been seen in 
the Hall resistance). A scenario resulting in such behavior can be constructed by noting 
that the liquid state will be the ground state only in a small region around 1/7; somewhat 
away from 1/7 the Wigner crystal state will become the ground state due to a cusp in 
energy at z/ = 1/7. Due to disorder induced spatial fiuctuations in u, the liquid state will 
therefore be the ground state only in a region of space where z/ is sufficiently close to 1/7. 
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Therefore, except for very small disorder, the WC is expected to percolate, but lakes of the 
incompressible liquid states should provide a reduction of the longitudinal resistance. At 
sufficiently small disorder, the regular FQHE behavior with well established plateaus may 
be expected. 

Finally, we come to the exciton dispersion at 2/5. The CP wave functions are more 
complicated here, so only 2 x 10^ the Monte Carlo have been performed for each energy, 
hence the larger statistical error in /Y^ in Fig. There are two strong minima at kl ~ 1.5 
and kl 0.75, consistent with our conjecture in Ref. [|l5l, based on the analogy to the IQHE, 
that the dispersion of the exciton of the n/{2pn + 1) FQHE state has n principal minima. 
We emphasize that even the biggest spherical system studied previously (with = 10 
particles [|12|) did not produce either the second minimum or the rise in the energy at small 
k. In general, larger systems are required for a reasonable description of the large-n FQHE 
states, since the relevant number of particles here is the number of composite fermions in 
the topmost CF-LL, which is (approximately) only Ne/n. A'^^ can similarly be computed 
for other FQHE states. We reiterate that no instability to a WC state is expected as one 
approaches u = 1/2 along the sequence n/{2n + 1), since the energy of the incompressible 
CF state remains substantially lower than that of the WC (Fig. 

Aside from the qualitatively new physics of instability, our study provides improved 
quantitative estimates for the excitation energies. The Table V quotes the thermodynamic 
estimate for the minimum energy required to create a CF exciton, A'^-^, for several FQHE 
states, along with the predictions using the single mode approximation. (Since both ap- 
proaches use the same ground state - the Laughlin wave function - the results indicate 
a lower energy for the CF exciton wave function compared to the SMA wave function.) 
The exciton energy is of experimental significance. The Raman scattering experiments are 
capable of directly measuring the energy of the exciton at the minima (or maxima) of the 
dispersion due to a peak in the density of states here |^ (although some disorder is required 
to provide a coupling to these modes, otherwise forbidden due to wave vector conservation). 
In addition, being the lowest energy excitation, A^j„ is a relevant parameter for a number 
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of other issues, such as the specific heat, phonon absorption p8|, and the finite temperature 
competition between the hquid and the WC phases |^ . 

As noted in the dispersion curves in Figs. become questionable at small kl, 
where a state with two excitons, each with energy A^j„, will have lower energy. There 



is some question as to whether the Raman peak associated with the kl ^ excitation of 
the 1/3 state is a single-exciton mode or a two-exciton (bound) state p2| , p9i , since the 
energy of the single exciton at /c/ ~ is close to 2A^j„ [^. At z/ = 1/5, on the other hand, 
the energy of the single-exciton mode at small kl is approximately 5A^j^, which should 
make it easier to distinguish it from a two-exciton state. 

In summary for this section, we have demonstrated that the CF scheme is now demon- 
strably capable of making detailed quantitative predictions, and, in particular, has sufficient 
accuracy to explain the collapse of the FQHE at small u, which manifests itself through a 
finite wave vector excitonic instability. 



IV. COMPOSITE FERMIONS WITHOUT HIGHER LANDAU LEVELS 

As discussed above, the complex system of strongly correlated liquid of electrons in the 
lowest LL can be understood simply in terms weakly interacting composite fermions moving 
in a reduced effective magnetic field and occupying, in general, several Landau levels of their 
own. Most of the essential features of experiments can be explained from this very simple 
picture, so it is no surprise that the microscopic formulations of composite fermions p|,^j2T[| 
have closely followed this physical picture. In particular, the wave function of composite 
fermions is obtained by taking the wave function of non-interacting electrons in the reduced 
magnetic field, occupying several LLs in general, attaching an even number of vortices to 
each electron to convert it into a composite fermion, and projecting the resulting wave 
function on to the lowest LL. Our confidence in this approach, referred to as the "standard 
approach" below, is further bolstered by the excellent quantitative agreement between these 
wave functions and the exact eigenstates. 
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Successful as this theory is, one cannot help but feel intrigued by the paradoxical feature 
that the analogy to higher LLs seems to play a central role in explaining the physics of 
electrons confined largely to the lowest LL. To be sure, there is nothing conceptually wrong 
with that. FQHE can and does occur at finite magnetic fields, when electrons are not fully 
confined to the lowest LL. And even in the limit of very high magnetic fields, when the 
electrons are strictly in their lowest LLs, all that one has done is travel somewhat outside of 
the allowed Hilbert space to a point where the relevant physical correlations are manifest, 
lock in this physics, and then continue back to the original Hibert space hoping on the way 
that it is not lost. In fact, analogy to higher LLs is one of the great strengths of the composite 
fermion theory since it facilitates a unified description of the IQHE and the FQHE and allows 
a back-of-the-envelope understanding of the key experimental facts without recourse to any 
microscopic theory whatsoever. 

Nonetheless, it would be desirable to understand composite fermions within the lowest 
LL. The challenge is not to write the CP wave functions in the lowest LL - the standard CP 
theory already accomplishes this - but to understand how to obtain within the lowest LL 
composite fermions along with their effective field, their effective LL structure, and effective 
kinetic energy, concepts that appear at the moment to be necessarily tied to the analogy 



with higher LLs. LLL composite-fermion wave functions have also been obtained in Ref. |^ 
from a different perspective, which does not capture the intuitive physics associated with 
composite fermions in an obvious manner. 

The LLL form of the CP wave functions that we described above sheds light on this 
issue. In order to emphasize this point, we will now show how the CP theory can be made 
plausible entirely within the lowest LL, i.e., develop what will be termed below a "LLL 
approach" for the CP theory, and answer the following questions: What is the structure of 
CP wave function in the lowest LL and how does it originate from the repulsive interelectron 
interactions? How does it produce an effective magnetic field? What is the meaning of the 
effective cyclotron energy and higher-LL-like physics in the lowest LL? In the end, we will 
show that the LLL approach is identical to the standard one. 
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For this purpose, it is more convenient to work in the disk geometry (i.e., the planar 
geometry with circular gauge), to which we will now switch. This will also serve as an 
extension of our new projection approach to the disk geometry, the consideration of which is 
useful in certain situations, e.g., for two-dimensional quantum dots in high magnetic fields. 



A. Basis states for a single electron 

In the circular gauge, the single particle eigenstates of an electron are given by 

Vn,m = iV„,„e-il^l^"L™(^) , (38) 

where z — x — iy denotes the position of the electron, the magnetic length is set equal to 
unity, the normalization coefficient 



/ n\ 

^"■"^^ V27r2-(n + m)! ' ^^^^ 
n = 0, 1, 2,... is the LL index, and the angular momentum is given by m = —n, —n + 1, ... 
In the lowest LL (n = 0), the single particles basis states have a particularly simple form, 
given by (omitting normalization factors) 

77o,m(^) = e-il^l^- . (40) 



B. Basis states for a single composite fermion 

Now, let us consider an (unsubscripted) electron at z in the presence of other elec- 
trons at Zj (for a total oi N -\- 1 electrons), which, at the moment, are treated simply as 
repulsive scatterers, and ask how the above basis states may be modified. Our aim is to 
write correlated basis states of the form 

e-il^l^'-Fl^,!^,}], (41) 

where the function F incorporates the effect of a repulsion between z and the other electrons 
Zj. The only way our electron can avoid the jth electron is through a factor like {z — Zj), 

28 



i.e., by creating a zero at Zj. (Factors of the type \z — Zj\°', which are nonanalytic in z, are 
not allowed due to the lowest LL restriction.) The first guess would be 

Fo[z.{z,}] = \{{z-z,) (42) 

j 

which has zeros of precisely located at the positions of other electrons, ensuring that 
z avoids all electrons. However, to deal with the general situation, we must find ways of 
constructing F[z^ {^j}] with different numbers of zeros of z. It might seem most straight- 
forward to create additional zeros, but this leaves us with too much freedom. We therefore 
resort to removing some of the zeros of Fq. First consider removing one of the zeros. The 
zero at Z}^ can be eliminated by writing {z — Zk)~^ — Zj). We symmetrize this function 
with respect to Zk and write 



Fi[zA^3]] 



. k 



Wiz - z,) . (43) 
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This function of z has — 1 zeros, which are not at the positions of the particles 
be seen by noting that one of the terms does not vanish under the substitution z = Zy 
However, since {2;^}] a/most vanishes when A^ — 1 of the terms vanish, it is 

expected that the zeros of this function are typically close to the particle positions. We have 
confirmed this numerically. Fig. ^ shows the positions of zeros of F\ (marked by circles) for 
a given configuration of Zj (crosses) for A^ = 200. There is a deficiency of a single zero near 
the origin, but away from this region the zeros are tightly bound to electrons. Thus, the 
removal of a zero shows up as a "defect" in the sea of electrons dressed with zeros. A study 
of systems of up to A^ = 1500 shows that the geometric average of the distance between 
an electron and the associated zero is given by ~ 1.43roA^~°'^'''^, where tq (= R/\/N for 
the droplet of radius R) is the typical interparticle separation, indicating that the zeros of 
Fi[z\ localize on to particles in the thermodynamic limit. Generalizing, n zeros are removed, 
creating n defects, by writing 



Fn[z,{Zj}] 
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n(^ - , (44) 



where the sum {ki, ki, ...kn} is over all distinct n-tuplets. Apart from an unimportant overall 
factor, Fn can also be expressed in the following equivalent forms: 



3 



lliz - z,) (45) 



F^[z,{zj}]^d-l[{z-zj), (46) 

3 

where d = d/dz and dj = d/dzj. 

The form of Fn[z, {%}] multiplying the single electron wave function satisfies, and in fact 
is completely determined by the following requirements: 

(i) It is restricted to the lowest LL. 

(ii) It has N — n zeros. 

(iii) It is a symmetric function of Zj. 

(iv) It is a homogeneous function of z and Zj, i.e., F„[az, {azj}] — a^~"'Fn[z, {zj}]. In 
other words, it is an eigenstate of the total angular momentum. 

(v) It is invariant under a translation of all coordinates by the same amount. 

The above scheme allows us to construct F with to zeros. Situations with more 
than N zeros can be handled straightforwardly by writing 

F,A^, {zj}] = - 'kT-'d^ Hi' - 'k) (47) 

3 k 

in which (p — 1)N zeros are bound exactly to electrons, and the other N — n zeros are the 
same as before. The zeros are more tightly bound to electrons in this function than another 
possible choice 

d^lUz-ZkY; 

k 

however, it is likely that both are valid. 

We define the basis functions of a composite fermion as (with the coordinates of other 
electrons treated as parameters): 

vZmi') ^ e-il^l%-+" - 'kT-'d- Hi' - 'k) , (48) 
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with m = — n, —n + 1, .... These are nothing but correlated basis functions of an electron, 
with some knowledge of repulsion between the electron under consideration and the other 
electrons explicitly built into them. They clearly are not orthogonal, but are in general 
linearly independent, so an orthogonal basis may be constructed. (Indeed, if we average 
over Zj, Vp,n,mi^) Vp,n',m'i^) orthogoual except for m = m'.) 

What is the energy of a composite fermion in state ripn^^{z)7 If we assume that the 
defects are independent, the energy of a composite fermion is proportional to the number 
of defects. We assign an energy of nEo to a composite fermion with n defects, where Ed 
is the (unspecified) energy of a single defect (with appropriate averaging over zj). This will 
be referred to as the independent-defect approximation (IDA). We will also see below that 
the simple act of counting the number of defects gives a considerable amount of insight into 
the energetics of the problem through the IDA with a minimal amount of work. Further, 
the Coulomb energies of the CF wave functions, which can be computed by Monte Carlo, 
will be shown to be excellent. 

An attractive feature of the single CF wave functions is that away from the defects the 
zeros are bound to electrons; i.e., there are no "wasted zeros". The no-wasted-zeros was 
an important property of the Laughlin wave functions for the ground states at LL fillings 
of u = l/(2p + 1) [^, which, however, could not be extended to the many-electron wave 
function at other fractions, even as a matter of principle. The present arguments shows how 
it can be restored quite generally, to an extent, at the level of the single CF wave function. 

Finally, defects at Ci, ■■■,Cn are created by replacing 9" by nr=i(^ ~ CD iii Eq. (^8]). 
The resulting wave function will be a linear superposition of rjp n.mi^)- ^^e explicit ex- 
amples below, we will only consider compact states composite fermions which will have no 
dependence on the positions of the defects after antizymmetrization; i.e., the replacement 
mentioned above will leave the full wave function unchanged. This is analogous to the wave 
function of an integer number fully occupied LLs of electrons: one may start with Wan- 
nier wave functions for electrons localized at (j, but the antisymmetrized many-body wave 
function is independent of Q. 
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C. Effective magnetic field 



One of the central features of the CF theory is the concept of effective field or the effective 
filling factor of composite fermions. This is a consequence of the fact that the zeros in the 
lowest LL are actually vortices, i.e., have phases associated with them. Let us take z in a 
closed loop enclosing a large area A and ask what is the phase associated with this loop. 
Aside from the Aharonov-Bohm phase, 27rAB/(f>o, there is a phase of —p27rpA due to the 
fact that z encloses ppA vortices in the loop, each contributing a phase of — 27r. (Here we 
neglect order unity corrections due to the presence of defects and an incomplete binding of 
the electrons and zeros. To see the minus sign, recall z — x — iy.) This is only half of the 
story though - it is important to remember that, in the final theory other electrons will also 
see p vortices at z, i.e., there is an additional contribution to the phase due to p vortices 
going around pA electrons, which is also equal to —p2'KpA. Equating the net phase to the 
Aharonov Bohm phase due to an effective magnetic field 5*-^^, 2t: AB'-^ ^ /4>o, gives 

B^^ = B- 2pp0o . (49) 

The filling factor of composite fermions, u'"^ = p4>o/B'-''^ , is related to the electron filling, 
u — p(l)o/B, as 

^ = ■ (50) 



D. Many-composite fermion states 

The basis states for many composite fermions are, as usual. Slater determinants of the 
single CF states: 
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CF 



(51) 



where the subscript of rjj denotes collectively the quantum numbers rij and mj of the com- 
posite fermion state (for a fixed p). The total angular momentum of this state is given 

by 



M = M^^ + pN{N - 1) 



where M*^^ = Y,f=i'n^j and the last term comes from Fpn- The total IDA energy is 



E = NdEd 



(52) 



(53) 



where No = H^=i % is the total number of defects. 



E. Testing the LLL wave functions 

The many-CF wave functions constructed above are trial wave functions, whose valid- 
ity remains to be confirmed. We again resort to studies of small systems for which exact 
solutions can be obtained numerically by an exact diagonalization of the Coulomb Hamilto- 
nian in the subspace of the lowest Landau level. We work in the planar geometry with free 
boundary conditions - confinement is achieved here simply by fixing the angular momentum. 
The exact ground state energy of seven electrons interacting via the Coulomb interaction is 
plotted in Fig. |T0| as a function of M. 

We begin by asking what insight can be gained simply by counting the total number 
of defects. It is clear that some defects will be necessary for angular momenta M^^ < 
^N{N — 1). It is straightforward to determine the minimum number of defects that N 
composite fermions must possess in order to produce a given total angular momentum of 
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M^^ = M — N{N — 1) (we take p = 1 here). The IDA ground state energy NdEd is also 
plotted in Fig. |T0|, with = 0.125e^/e/, chosen empirically to get the best fit. There 
is a remarkable similarity between the shapes of the two curves and the structure of the 
cusps, providing strong support for the qualitative validity of the CF picture and our IDA 
association of the energy with the number of defects. 

We next carry out a more quantitative test of the microscopic structure of the CF wave 
functions. The number of IDA ground states at M'"^ is tremendously small compared to 
L, the total number of electron states at a given M. The most dramatic reduction of the 
low energy Hilbert space due to the formation of composite fermions occurs at those values 
of and M for which there is a unique IDA ground state. The CF ground state here has 
a compact pyramid-like structure: the Nn composite fermions with n defects occupy the 
smallest angular momenta (i.e., m = —n, —n + 1, ... —n + N — 1), and Nq > Ni > N2.... The 
wave function of such a compact many-CF state is a single Slater determinant, denoted 
by [A^o, A^i, ...]. We have computed the Coulomb energies of a number of compact CF 
ground states using Monte Carlo techniques. As seen in Fig. [10, the results are virtually 



indistinguishable from the exact ground state energies. Table V shows the CF energies for a 
number of compact states for up to seven particles. The deviations from the exact energies 
typically appear in the fourth significant digit, and are small compared to the differences 
between energies at different M. These results provide a substantial, non-trivial test of 
the CF theory for two reasons. First, the wave functions of compact CF states contain no 
adjustable parameters; they are completely determined by symmetry within the space of 
the CF wave functions. Second, the size of the many- electron Hilbert space, L, (which is 
the dimension of the matrix diagonalized to obtain the exact ground state) is in general 
quite large, as shown in Table V. (This was incidentally also true for the examples discussed 
earlier in the context of the spherical geometry.) The only exception is at M = A^(A^ — l)/2, 
for which there is one one state in the lowest LL, the u = 1 state. This corresponds to the 
[1, 1, ...] of composite fermions. In this case, the LLL CF wave function is exact for trivial 
reasons, as also confirmed explicitly in our calculations. Note also that when there is more 
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than one compact state possible for a given M (e.g., for = 6 and M = 13 in Table I), the 
one with the smallest energy ought to be identified with the true ground state. 



F. Equivalence between the standard and LLL approaches 

The discussion above makes no reference to higher LLs, yet obtains the effective magnetic 
field of the composite fermions. In this section, we show that it is in fact identical to the 
standard approach. In the latter, the unprojected wave function of the many CF state, 
(^CF-up^ is related to that of the non-interacting electrons, "I*, as 



^CF-up ^ -Q , 



|2P $ 



(54) 



j<k 



(We note here that the gaussian factor in $ is taken at the external magnetic field. Alterna- 
tively, $ can be written at the effective magnetic field, but then a gaussian factor must be 
supplied along with the Jastrow factor, which, for bulk FQHE states, implies multiplication 
by where $i is the wave function of one filled LL 0.) 

For illustration, consider only states for which $ is a single Slater determinant, made up 
of single electron states of Eq. (0). Then, 



j<k 



^2(ri) r/2(r2) 



(55) 



Here the subscript of rjj denotes collectively rtj and mj, the LL index and the angular 
momentum, respectively, of the electron. 

The total angular momentum of $'^'^~"p is given by 



M = M^^ +pN{N -I) 



(56) 



where M^^ = J2i 

rrij is the angular momentum of fermions in state $. In the mean- field 
approximation (MFA) of the standard picture, the composite fermions are treated as non- 
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interacting, so their total energy is equal to the kinetic energy of electrons in the state 
with the cyclotron energy replaced by an effective cyclotron energy of composite fermions; 



I.e., 



E = nj]huj 



CF 



(57) 



The equivalence between the the LLL and the standard schemes follows if we make the 
following identifications: 



number of defects = LL index 



(58) 



anticipated by our notation, and 



energy of a defect (Ed) = effective cyclotron energy (hu 



CF\ 



(59) 



The IDA of the LLL approach is then identical to the MFA of the standard scheme. In 



particular, the curve of the defect energy of composite fermions in Fig. |T0| is identical to 
that considered earlier in [Hl| for the effective kinetic energy of noninteracting fermions at 

We now establish the equivalence between the two approaches at the level of the micro- 
scopic wave functions themselves. Noting as before 



j<k j 



(60) 



where the prime denotes the condition j ^ k, the wave function can be written as 



(^CF-up 



CF-up, N CF-up, ^ 



CF—upr \ CF—upr \ 



(61) 



where 
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riZrr^n^^) = VuA^^) n(^. - (62) 



is interpreted as the unprojected wave function of the jth composite fermion |^5|. The CF 
wave function in the lowest electronic LL is obtained by a LLL projection of <|)'^^~"p, which 
can be accomplished by projecting rf'^-'^v lowest LL. 

The crux of the argument is that a LLL projection of f]pn^^ produces precisely (apart 
from normalization factors) rip^n,m ^'i- (El)' 



i k 

This shows that of the standard approach is identical to of the LLL approach. 

To prove Eq. (|63|), we first write 

where /cq = 'max(0,— m). Our goal is to evaluate the projection of the product 
rin,m{z, z*)Y[k{z — Zk)- As showu in Refs. [^,0 the projection is obtained by first bring- 
ing all 2*'s to the left of 2;'s in the polynomial part multiplying the gaussian factor and then 
making the replacement z* 2d , while remembering that the derivatives do not act on the 
gaussian factor. Le., 

^^n,m(r) n(^ - Zk) = VnA'^) Y[{z " Zk) , (65) 

k k 

with the operator rj given by 



k=0 



Substituting 



and using 



Qk^k+m _ kl^fk + m\ 



n k n n 

E E = E E 

k=ko a=ko a=ko k=a 
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the sum over k is seen to be proportional to (1 — 1)" " which vanishes unless a = n. Thus, 
only the term with k = a = n survives, giving for the operator r/: 

r7,,^(r) = iV„,^izl)!e-il^l'z™+«5" . (67) 

Substitution into Eq. ( |63D reproduces Eq. p8|), completing the proof of equivalence of the 
two formulations of the CP theory. 

The summary of this section is that we have shown that the origin of composite fermions, 
their effective magnetic field and LL structure, as well as their wave functions can all be 
obtained rather straightforwardly within the lowest LL. This complementary description 
gives additional insight into the physics of composite fermions, the formal structure of their 
wave functions, and the role of the LLL projection operator in the standard approach. 

V. CONCLUSION 

We have developed an explicit LLL representation of the CP wave functions in which 
the LLs of composite fermions are filled in the standard manner using the single particle 
basis states of composite fermions, in complete analogy to the problem of non-interacting 
electrons. The most important consequence is that it is now possible to study large CP 
systems, containing up to 50 composite fermions. Given how much insight has been gained 
from the study of systems with ten or fewer particles, this represents a significant setp toward 
a better quantitative understanding of the PQHE. In particular, certain interesting regions 
of filling factor have become accessible for the first time to quantitative investigations. We 
have obtained estimates for the ground state energies and gaps (both to charged and neutral 
excitations) for a number of PQHE states. The charge gaps predicted by the CP theory 
are consistent with the general behavior observed experimentally and that anticipated by 
a mean-field theory, although more work will be needed to include the effects of non-zero 
transverse width and LL mixing before the theoretical gaps may be compared directly to the 
experimental ones. An important fact discovered in the course of this work is that the CP 
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theory is capable of explaining the absence of FQHE at small filling factors (z/ < 1/9), where 
the FQHE is explicitly shown to be unstable to a spontaneous generation of CF-excitons. 

The new LLL representation of the CF wave functions also helps resolve the longstanding 
question: do we necessarily need to use the analogy to higher LLs to explain the physics 
of electrons confined to the lowest LL? The answer is in the negative. We have shown that 
Laughlin-like arguments can be used to obtain CF wave functions within the lowest LL of 
electrons. Concepts like the CF-LL index, CF cyclotron energy, CF filling factor, etc., which 
seem to be inextricably tied to the analogy with higher LLs, appear naturally also in our 
LLL approach. In particular, the CF LL index appears as the number of defects in the 
single-CF wave function and the CF cyclotron energy as the energy of a single defect. 

This work was supported in part by the John Simon Guggenheim Foundation, the Na- 
tional Science Foundation under grant no. DMR9615005, and by the MRSEC Program of 
the National Science Foundation under award number DMR 94-00334. We thank Tetsuo 
Kawamura, Daniil Kaplan, and K. Park for discussions and help with numerical work. 



n\f 2q + n \ / 5 \ ' .q+a/+s / ^ ' " ' 



^Q+M+s ^Q~M+n-s 



VI. APPENDIX 

Here we prove the general correctness of the formula 

VY Y rn I = Ke^^'^^ 

I ^ q,n,m^ q' ,(),m' -ivo 
n 

' \sj\q + n — m — sj ^du J V^"^/ 
where M = m + m' and Q = q + q'. (Here, K represents factors that will not be relevant for 
the arguments below.) The validity of this equation was shown for |M| < Q in Section J|. 
In order for it to be valid generally, the right hand side must vanish identically for \M\ > Q, 
since no states exist in the lowest LL with \M\ > Q. This is what we now show. 
First consider the case M < —Q. This would be possible when 

m = —q — a, a = 1, ...n , 

m' = -q' + (3, /5 = 0,1,... 
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Then 



M = m + m'=-Q-a + /3, 



and the condition M < —Q imphes that a — (3 > 1. In fact a — (3 — 1, 2...n. In this range, 
the binomial coefficients impose the following constraints on s: 



As a result, the sum over s goes over 

s — a,a + 1, ...n. 

Now, consider the derivative (^^j ^Q+m+s r^^^e exponent of u can be written as 

(5 + M + s = s-(a-/5). 

It is always less than s, since a — j3 >1, and also non-negative, since the smallest value of 
s is q; (> 1 by definition). Therefore, the derivative with respect to u vanishes: 



s < n 



s<q + n — m — 2q + n + a 



s>0 



s > —m — q = a 




It can similarly be shown that for M > Q, 




V 



,Q-M+n-s _ g 
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TABLES 

TABLE I. The ground state energies per particle for several systems at = 2/5 and 3/7 are 
obtained from the CF wave functions (third column) and exact diagonalization (fourth column). 
The fifth and sixth columns contain the energies of the charged excitation, in which one composite 
fermion from the south pole in the highest occupied CF-LL is removed and placed on the north pole 
in the lowest unoccupied CF-LL. The seventh and eighth columns give the CF gaps and the exact 
gaps. The energies are in units of e^/eZ, where I is the magnetic length and e is the background 
dielectric constant, and include interaction with the uniform positively charged background. The 
statistical uncertainty in the last digit (s) of the CF energy is shown in brackets. The exact results 
for A = 10 and 12 are taken from S. He, S.H. Simon and B.L Halperin, Phys. Rev. B 50, 1823 
(1994). 
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TABLE II. The CF energies for ground states at several filling factors as a function of N. The 
charge gaps are also shown in each case; the corrected gap is obtained by subtracting the interaction 
energy of two point-like particles of appropriate charges at the opposite poles (see text) . The Monte 
Carlo uncertainty for the gaps is the same before and after correction. 
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TABLE III. The thermodynamic estimates for the ground state energy per particle and the 
charge gap for a number of FQHE states. 
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TABLE IV. The minimum gap to neutral excitation, A^^^, for several filling factors. The SMA 

prediction is also given for comparison. 
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TABLE V. Comparison between the energy of the CF wave function and the exact 

Coulomb energy for the ground states at several angular momenta M in the disk geometry. The 
energies are in units of e^/el. No positive background has been assumed; confinement is due to 
fixing the total angular momentum. The energies of composite fermions in the compact states 
[A^o, Ni,...] (notation explained in text) have been evaluated by Monte Carlo, with the statistical 
uncertainty in the last few digits shown in brackets. The exact energies have been obtained by a 
numerical diagonalization of the Coulomb Hamiltonian in L dimensional Hilbert space. 
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FIG. 1. Ground state energy of various FQHE states as a function of the number of particles, A^. 
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FIG. 2. The energy of the hquid FQHE states and the correlated WC state. The WC energy 
is taken from Only the quantum correction to the energy of the classical ground state (a 

two-dimensional Wigner crystal with triangular symmetry) is shown, where E^-i = — 0.782133-^ 
e'^ /el. The open diamond on the right vertical axis is our estimate of the CF ground state energy at 
u = 1/2, — Eel = 0.0877(2)e^/e/) obtained by extrapolating the energies of the n/(2n -|- 1) 

CF states. Also note that the energy of the CF liquid is shown only at the special n/{2pn + 1) 
filling factors; the full curve will have cusps at these points. 
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FIG. 3. N dependence of the "corrected" gap for 2/5 and 4/9. 
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FIG. 4. Thermodynamic values of the gaps plotted as a function of l/(2n + 1). The results for 
the 1/3 gap are taken from Ref. The straight line is obtained by a chi-square fit, given by 

Eg = 0.000(5) + 0.316(21)(2n + 
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FIG. 5. Exciton dispersions for FQHE states at v = 1/5 for systems with several sizes, with 
N shown on the figures. All curves assume zero thickness. The error bars are shown explicitly 
whenever they exceed the sizes of the symbols. 
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FIG. 6. Same as in ^ for i/ = 1/7. 
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FIG. 7. Same as in ^ for i/ = 1/9. 
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FIG. 8. Exciton dispersions for i/ = 2/5 for N = 30. The dashed hne is a guide to the eye. 
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FIG. 9. Positions of zeros (dots) of {zj}] for a given configuration of zj (crosses) for 
A'^ = 200. The other electrons {zj) are taken to be in a reasonably probable configuration in the 
presence of a repulsive interaction. A defect is seen near the origin. 
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FIG. 10. The empty circles show the exact ground state energy of seven interacting electrons 
confined to the lowest LL in the angular momentum M sector obtained by an exact diagonalization 
of the Coulomb Hamiltonian. The filled triangles give the independent-defect-approximation for 
the energy of composite fermions in M^^ state, NdEd, where Nj:) is the minimum number of 
defects at and Ed is the energy of a single defect, chosen to be Ed = 0.125e^/e/ to obtain 

the best fit (also, an arbitrary overall constant has been added to the defect energy). The filled 
dots show the Coulomb energies, evaluated by Monte Carlo. The configurations of the compact CF 
states are indicated on the figure in sequential order. The wave functions of these states contain 
no adjustable parameters. The lines are a guide to the eye 
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